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A theory for the effects of nonmagnetic disorder on the magnetic pair breaking rate a induced 
by paramagnetic impurities in quasi two-dimensional superconductors is presented. Within the 
framework of a strong-coupling theory for disordered superconductors, we find that the disorder 
dependence of a is determined by the disorder enhancements of both the electron-phonon coupling 
and the spin-flip scattering rate. These two effects have a tendency to cancel each other. With 
parameter values appropriate for Pbo.gBio.i, we find a pair breaking rate that is very weakly depen- 
dent on disorder for sheet resistances < Rn J$ 2.5 kS7, in agreement with a recent experiment by 
Chervenak and Valles. 
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I. INTRODUCTION 

The physics of the observed T c suppression in super- 
conductors that contain nonmagnetic disorder has been 
the subject of much debate in recent years. Let us fo- 
cus on homogeneously disordered thin superconducting 
films, with the disorder parametrized by the normal-state 
sheet resistance Rn, and on the BCS-like quasi tran- 
sition that is well pronounced in these films although 
the true superconducting transition is of Kosterlitz- 
Thouless natureE In these systems, the BCS transi- 
tion temperature T c , defined as the mid-point of the 
resistive transition, is observed-to decrease monotoni- 
cally with increasing disorder.Em A complete quantita- 
tive understanding of this effect within a microscopic 
strong-coupling theory has proven difficult, although 
the first perturbatiiae calculations within a phenomeno- 
logical BCS modelO were rather promising. Qualita- 
tively, disorder-induced changes in the electron-phonon 
coupling!} in the Coulomb repulsion between the con- 
stituents of theXlooper pairsja and in the normal-state 
density of statescl have all been identified to be important. 
The difficulty lies in the fact that some of these effects 
tend to suppress T c while others tend to enhance it, and 
T c depends exponentially on all of them so that subtle 
balancing effects result. Also, the number of parameters 
that acquire a disorder dependence is quite large, and fits 
of theoretical results to experimental data are therefore 
not necessarily very conclusive. Indeed, theories that are 
structurally quite different, and mutually inconsistent, 
have been shown to fit the same T c data equally well.Lrla 

In this situation it is obvious that one should study the 
disorder dependence of quantities other than the transi- 
tion temperature in an attempt to discriminate between 
various theories, and to obtain independent information 
about the disorder dependence of the various parameters 



that determine T c . One possibility-is to measure the in- 
elastic lifetime of quaaparticlesEllj The experiment by 
Pyun and Lembergeul3 on amorphous InO has been an- 
alyzed by the present authorsEU in the framework of-a 
strong-coupling theory for disordered superconductors ,E3 
and quantitative agreement between theory and experi- 
ment has been achieved. Another possibility is to study 
the influence of nonmagnetic disorder on the pairbreak- 
ing induced by a small amount of magnetic impurities in 
addition to the nonmagnetic ones. This has the advan- 
tage that the pair breaking parameter is easier to measure 
than inelastic lifetimes, and that it can be measured si- 
multaneously with the T c suppression in a series of films 
where the nonmagnetic disorder is varied in situ by con- 
trolling the film thickness. 

Such an experiment—has recently been performed by 
Chervenak and Valles ,E3 who studied quench condensed 
ultrathin films of Pbo.gBio.i of varying degrees of disor- 
der (150 < Rn < 2.2 kf2, leading to a T c between 6 
K and 2.35 K). Of each sample, one half was doped with 
Gd, while the other half was left undoped. Gd acts as a 
paramagnetic impurity and leads to pair-breaking and a 
reduction of T c . The transition temperature as a func- 
tion of Gd concentration was then studied as a function 
of film thickness, which is correlated with R D . The re- 
markable result was that the pair-breaking parameter a 
is only mildly dependent on disorder for films with nor- 
mal state sheet resistances Rn ranging from 0.15 kSl to 
2.2 kfi. The implication seems to be that the effects 
of disorder that lead to lower transition temperatures 
do not manifest themselves in the spin-flip pair break- 
ing rate. An attempt to understand this behavior by a 
phenomenological modification of the Abrikosov-Gorkov 
result for a, using the renormalization of the density of 
states inherent in Refs. HJlJ and ||, failed.E3 This poses 
the important question whether the success of these the- 
ories in describing the T c suppression and the disorder 
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and temperature dependence of the inelastic lifetime was 
fortuitous, and whether they are lacking some important 
physical ingredient that manifests itself in the pairbreak- 
ing rate. 

It is the purpose of the present paper to analyze these 
questions. What we will find-is that one needs to take the 
strong-coupling correctionst-3 to the Abrikosov-Gorkov 
expression into account before one generalizes to the dis- 
ordered case in order to obtain the correct structure of 
the theouy. Once this is done, we find that our previous 
theorycaO accounts very well for the observed effect. 



II. FORMALISM 

Our starting point is aux theory for the suppression of 
the superconducting T c x2 and the enhancement of—the 
inelastic scattering ratejLil by nonmagnetic disorder S3 In 
this section we recall the most salient features of this 
theory. First one uses an exact eigenstate formalism 
to derive generalized Eliashberg equations for the nor- 
mal Green function G(e, iui), and the anomalous Green 
function F(e, iui). Since the wavenumber is not a good 
quantum number in the presence of static impurities, G 
and F are functions of energy and frequency, rather than 
wavevector and frequency as in Eliashberg theory. The 
Green functions are expressed, as usual, in tcrms^of an 
anomalous self-energy, W(e, iui), and a normal onc.E3 The 
latter is split into a piece iuiZ(e, iui) that is an odd func- 
tion of frequency, and a piece Y(e, iui) that is even in ui. 
In terms of these quantities, G and F read, 
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-W(e,iu>) 



[iuiZ{e,iui)} 2 - [e + Y{e,iui)} 2 



(la) 



(lb) 



In clean systems, the normal self-energy piece Y is a con- 
stant that just shifts the chemical potential and can be 
omitted. In the presence of disorder, pipapyer, Y has 
been found to be of crucial importance,t30ll3 and to re- 
flect ther-ahysics of the Coulomb anomaly in the density 
of stateal3 in the context of superconductivity. The gen- 
eralized Eliashberg equations then take the form of inte- 
gral equations in both energy and frequency for the three 
functions Z, W, and Y. A solution of these equations has 
been obtained by means of various approximations. In 
particular, the frequency dependence of Y was found to 
be weak and could be omitted, and its energy dependence 
was approximated by the first term in a Taylor expansion 
about a characteristic energy ui, 



Y{e,iu>) ~ (e-Q)Y' 



(2) 



with Y' r^f|— ,dY/de\ e —Q, and u> an average phonon 
frequency.c2Ea With some further approximations, the 



energy integrations could then be performed, and the 
theory be cast in the same form as Eliashberg theory. 
A two-square well approximation then leads to a T c for- 
mula that has the structura-of a generalized McMillan or 
Allen-Dynes formula, vizJl3iI3 
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exp 



-1.04(1 + X + Y') 



A-/2*[1 + 0.62A/(1 + F') 



(3) 



Here Y' is the normal self-energy piece mentioned above, 
and A and fx* are disorder dependent generalizations 
of the electron-phonon coupling constant A and the 
Coulomb pseudotential /i*, respectively, in Eliashberg 
theory. Explicit expressions for all three of these quanti- 
ties have been given in Ref. [l2|, and will be evaluated for 
the case of thin films below. 

In the presence-of magnetic impurities, T c is reduced 
by pair breaking^ Abrikosov-Gorkov theor*-,has been 
modified to allow for strong-coupling effects,E-3 with the 
only resulting change being a factor of 1/Z in the pair 
breaking parameter. The result is 

- ln(T c /T c0 ) = ^{a/2nT c + 1/2) - ^(1/2) , (4) 

with T c o the value of T c in the absence of the magnetic 
impurities, tp the di-gamma function, and a the pair 
breaking parameter a = (1/Z)(1/t s ), with l/r s the spin- 
flip scattering rate. In the case of clean superconductors, 
Z = 1 + A. In the presence of nonmagnetic disorder, and 
with the same approximations that lead to the T c for- 
mula given by Eq. (^), it is straightforward to repeat the 
calculation of Allen and MitrovicE-3 within the framework 
of the theory of Ref. |l2|. The result is again Eq. (Q), but 
with a replaced by a disorder dependent a which in turn 
is related to disorder dependent parameters, 



a 



1/r. 
1 + A 



(5) 



Here A is the same quantity as in Eq. (||) , and l/f s is the 
disorder dependent spin flip-scattering rate. Throughout 
this paper, we choose units such that h = ks = 1- In 
the next section, we derive an explicit form for a in a 
disordered thin superconducting film. 



III. DISORDER DEPENDENCE OF a 

A. Electron-phonon coupling A 

The electron-phonon coupling strength A is defined as 
an integral over the Eliashberg function a 2 F, 



(6) 



a 2 F, and hence A, are disewier dependent due to effects 
first discussed by PippardE3 in the context of ultrasonic 
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attenuation, and by Schmidcl for the electron-phonon in- 
elastic lifetime. The main physical point is that disorder 
decreases the coupling of the electrons to longitudinal 
phonons due to collision drag, but increases the coupling 
to transverse phonons due to the breakdown of momen- 
tum conservation. For realistic parameter values the lat- 
ter effect is stronger than t he-former, leading to an overall 
increase of A with disorder .E3 For Debye phonons in 3 — d 
systems, A has been calculated in Ref. 12. Repeating that 
calculation in d = 2 is straightforward, and very similar 
to the corresponding calculation-of the electron-phonon 
inelastic lifctimc.cil The result isE3 



A = 21 — £ § f b (ul/c b ) 

b=L,T b 



(7) 



for a system with mean free path /. Here cl,t are the 
longitudinal and transverse speeds of sound, respectively, 
Lj£> is the Debye frequency, the dimensionless constant 
db = kp/lGirpiCb with pi the ion density and kp the 
Fcjgii wave number, and the functions fp,L are given 
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(8b) 
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Substituting Eqs. (^) into Eq. (^) we obtain the disorder 
dependence of A, 



A = 



XRuEpcp 



ttujdVf 
RuEfCl 



+ 2-±F T 



TTLUp,Vp 



,R a EpCT 
where we have defined two functions, 



(9a) 



F L (x) = Vl + x 2 - 1 - ln[(VT 
1 - y/T+a? 



F T {x)-- 
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l)/2] , 
* + l)/2] 



(9b) 



2(1 

with A = Aiop)dhl'Kmc 1 L the electron-phonon coupling in 
a clean 2 — d system. The dimensionless resistance R a = 
Rne 2 /Ti ps Rn/4.1 kfl is a measure of the disorder in the 
material. 

As in three dimensions, the size of the disorder renor- 
malization of A depends on the ratio of the longitudinal 
to the transverse speed of sound. This is a result of the 
abovementioned competition between an increase in the 
coupling between electrons and transverse phonons and 
a decrease of the coupling to longitudinal phonons. Since 
the transverse speed of sound is invariably smaller than 
the longitudinal one, A increases with increasing disor- 
der. This effect tends to reduce the pair breaking rate, 
Eq. (5). However, we also have to calculate the disor- 
der dependence of the spin-flip scattering rate in order 
to obtain the disorder dependence of a. 



B. Spin-flip scattering rate l/r s 

The interaction between the electron spin and an im- 
purity spin S(r) at site r is described by a Hamiltonian, 



Ha 



k,k',/i,i/ 



Jk.k' S(k- k') ■ (typ&nuCiw 



(10) 



Here c' and c are fermion operators, a = (a x , a y ,a z ) de- 
notes the Pauli matrices, and Jk,k' denotes the electron- 
magnetic impurity exchange integral. We now calculate 
the electron self energy contribution, E, due to this inter- 
action in Born approximation. It is most convenient to 
do this in an exact eigenstate representation, in analogy 
to the calculation of the Coulomb self energy in Ref. |2^. 
The calculation is straightforward, and we obtain 



S(e, iuS) 



J ^G{e'MY,Vs{<l)F a {<l,e-e') .(11) 



Here, Np is the free electron density of states per spin at 
the Fermi level. We only retain the s— wave component of 
the interaction so that Vg(q) = npS(S+ I) J 2 , where np 
is the concentration of paramagetic impuritiesp-and J is a 
measure of the exchange interaction strengthBH G(e, iio) 
is the normal Green function in the superconductor and is 
given byEq. @. Finally, F s is the spin density analogue 
of the density-density correlation function denoted by F 
in Ref. |2^. If we work to lowest order in the electron- 
impurity spin interaction, and neglect Coulomb and finite 
temperature effects in S, then F and F s are identical. 

We obtain the spin- flip scattering rate l/r s from the 
self energy S by analytically continuing to real frequen- 
cies, iu> — > lo + iO, and going 'on-shell', i.e. putting 
e = lu. For our purposes, we are interested in the in- 
fluence of spin-flip scattering on the superconducting T c . 
The physics that determines the latter is dominated by 
processes on a frequency scale of w, a typical phonon 
frequency. For the same reason for which we take the 
parameter Y' in Eq. (|^) at the frequency to we therefore 
define 



1/t s = — 2ImS(a;, iu — > lu + iO) 



(12) 



In a clean system, the spin-density correlation func- 
tion F s (q,uj) is frequency independent. In AfX case we 
recover from Eq. (^2|) the well known resultEjEJ 



n P S{S+ l)J 2 4N F 



(13) 



In a disordered system, F s is diffusive]!] and in the Green 
function G we must keep the self energy piece Y' as dis- 
cussed above. We thus obtain 
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where 



2n P S(S+ I) J 2 
N F \l + Y'] 



q,uj 



A 



1 + Y' 



(14a) 
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F s (q,uj) = 5(g)- 



Dq 2 



{Dq 2 ) 2 



(14b) 



with g(q) the Lindhard function, which for simplicity 
we replace by NpQ(2kp — q). Here D denotes the nor- 
mal phase spin density diffusion coefficient, which in the 
nonintcracting electron approximation coincides with the 
mass or charge diffusion coefficient, so D = n/ mRu- Per- 
forming the wavenumber integral in Eq. ( 14a ) we finally 
obtain 




i+i4^ 



with 



to 



1 + Y' 



(15a) 



(15b) 



and 1/t s given by Eq. (|13J). l/f s depends on disorder 
both explicitly, and implicitly through Y'. Our final task 
is therefore to calculate the dependence of Y' on R n ■ 



C. Normal Self Energy Piece Y' 

In order to calculate Y' we again have to repeat the cal- 
culations of Ref. [[2] in d = 2. Both the electron-electron 
and the electron-phonon contributions to the self energy 
contribute to the self energy piece Y. Performing a Tay- 
lor series expansion in energy around e = wof Eq. (2.12) 
of Ref. n3, we obtain 



Y'(Q) = 6U%(Q)+4 / — 6a 2 F H (Q,v) 



dv 



(16) 



5Uq(lu), which describes the Coulomb contribution to 
Y', is taken from Ref. H, 



Dq 2 



Vc(q) 



5(q) : 



{Dq 2 ) 2 +U) 2 



(17) 



H 9k,k'(q)5p, P '(q)^c(k-p) 



k,p k',p' 

with the statically screened Coulomb potential, 

1 K 



Vc(q) 



2N F K + q 



K = 4ne z N F 



(18) 



Using the prescription for performing momentum sums 
of this type as described in Ref. [12], the integrals can be 
done and yield 
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8irE F k 2 ' 8irE F ' k 
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with the functions 
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The Coulomb pseudopotcntial /1 in d — 2 is given by 

K 



(20) 
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2k F 
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(21) 



As discussed in Ref. [12], the phonon contribution to 
Y' , which is given by the second term on the right-hand 
side of Eq. (|l6|), is related to a stress-stress correlation 
function and can be calculated in a similar manner as 
SUq . One obtains 



4 y ^Sa 2 F H (u,,v) = 

.A C L E F . _! 

Attn tt^ sm 

Finally, both contributions can be collected to give 
Y'(lu) = R 



(22) 
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(23) 



where G and H have the same arguments as in Eq. 
(|l9|). The enhancement of Y' with increasing disorder 
is due to the opening of a correlation gap in the (normal 
state),-deHfiity of states, and contributes to the decrease 

of TcMta 



IV. FINAL RESULT AND DISCUSSION 

We are now in position to collect our results and thus 
obtain the disorder dependence of the magnetic pair 
breaking rate. Substituting Eqs. @, ([l5|), and ( p3| ) into 
Eq. (H) yields our final result, 



a 
a 



1 + X 
1 + X 



1 R c 
1+Y' 8^ 



In 



8nE f 



□ or 



(24) 



The disorder renormalizations of the pair breaking rate 
appear both in the numerator and denominator and 
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therefore the rate may either increase or decrease with 
increasing disorder depending upon the material param- 
eters X,cl,ct,vf, Ef,ojd, and [i. 

We now address-the experiment on Pbo.gBio.i by Cher- 
venak and VallesO. To estimate the parameters enter- 
ing into Eq. (24), let us first consider the parameter 
values for bulk PbBi, as far as available, given in Ref. 
|(| Thereby we have ui = 56K, andu D = 108K. The 
Bohm-Staver relation gives 



E F c L 



2(2/ Z) 



1/3 



(25) 



For clean bulk Pb one has E F = 1.1 X 10 5 K, c L =2050 
m/s, ct — 710 m/s, and Z = 4. However, we do not ex- 
pect the actual parameter values to correspond to those 
for either bulk Pb or bulk Pbo.gBio.i. First of all, the 
material in question is a thin film, and moreover the sub- 
strate is expected to modify its properties, in particular 
the acoustic ones. Evidence for this is provided by the 
fact that the parameter values quoted above do not give 
the correct value of the clean limit T c as measured in Ref. 
|l3| . It is therefore likely that the substrate on which the 
thin layer of Pb is deposited strongly affects the phonon 
spectra of the film, altering both A and the ratio cl /or 
compared to bulk Pb. Accordingly, we choose a value for 
the bare A = 1.12 and fi — 0.1 which (in the absence of 
disorder) reproduces the highest T c as measured in Ref. 
|i~3| . Lastly, since cl/ct is not known even for Pbo.gBio.i, 
we let cl/cl be determined by a fit to the data. 

These parameters provide the curves shown in Figure 
1. The solid curves are the results for the disorder depen- 
dence of the normalized pair breaking rate as given by 
Eq. (5). The points represent the data taken from Ref. 
|i"3| for two different runs. The decrease at small Rn is due 
to the fact that for small disorder, A and Y' grow more 
rapidly than l/f s . The normalized rate goes through a 
shallow minimum at roughly Rn ~ 0.3 at which point the 
disorder renormalizations of Y', A, and l/f s are balanced 
and offset each other. With further increasing disorder 
the enhancement of l/f s dominates, and leads to a slowly 
increasing pair breaking rate. 

We remark that the point at which the minimum oc- 
curs depends sensitively on the ratio of the longitudinal 
and transverse speeds of sound. To obtain the solid lines 
in Fig. la, cl/ct — 1.9 was used while 2.1 was used for 
Fig. lb. These values lie between the value for bulk Pb 
(2.88) and the substrate (similar to pyrex, 1.72) used in 
Ref. |l3| and thus does not seem unreasonable. Larger 
values of cl / Ct yield a more drastic reduction of the rate 
for small disorder and the region of increasing a occurs at 
larger values of Rn . This sensitivity of the overall shape 
of the curve to the material parameters may be reflected 
in the slightly different results obtained for the two ex- 
perimental runs in Ref. [fl] as shown in Fig. 1 . It would 
therefore be very interesting to repeat the experiments 
using substrates with different acoustic properties. 

In summary, we have presented a theory for the param- 
agnetic pair breaking rate in disordered superconducting 



films and have shown that the disorder dependence of 
the rate depends delicately on the disorder renormaliza- 
tions of Y' , A, and l/f s . As a result, the rate can either 
increase or decrease with disorder, depending upon ma- 
terial parameters, and in general it is not a monotonic 
function of disorder. Our conclusion is that the disor- 
der dependence of the rate as observed by Chervenak 
and Valles in Pbo.gBio.i films can be quantitatively un- 
derstood via an application of the microscopic theories 
developed in Refs. |ll] and |i~2| . 
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FIG. 1. Plot of the pair breaking rate as a function of dis- 
order for two different Pbo.gBio.i films. The symbols are 
data from Ref. [l3|, the solid line is obtained from Eq. (|B|). 
For Fig. la(b) we have chosen cl/ct = 1.9 (2.1), and 
a = I / t s {1 + A) = 90.5K (70K). All other parameters are 
the same for both figures and are given in the text. 
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